% 创建时间：2025-06-20 23:18:18
% 作者：兪馳
clc;
clear;

A=[1 2 3;
    4 5 6;
    7 8 0]; % 输入已知矩阵
C=-[1+1i, 3+3i, 12+10i; 
    2+5i, 6, 11+6i; 
    5+2i, 11+1i,2+12i]; % 注意下面的转置
A0=sym(kron(eye(3),A)+kron(A,eye(3))); % 构造系数矩阵
c=C(:); % 系数矩阵与向量
x0=-inv(A0)*c; % 求解析解
x=reshape(x0,3,3) % 恢复解矩阵
norm(A*x + x*A.' + C) % 验证结果